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THE GEOCENTRIC ORIENTATION VECTOR 
FROM LIMITED ASTRO- GEODE TIC DATA 

R. S. Mather* 

Geodynamics Branch 


SUMMARY 

The establishment of a world geodetic system from gravity 
data available at the present time, can be based on the corrections 
necessary to convert the separation vectors, as determined from 
astro-geodetic information, to equivalent values which could be 
directly compared with those determined from an imperfect know- 
ledge of the Earth’s gravitational field. 

Test computations using the comprehensive sets of both astro- 
geodetic and gravimetric data available for the region covered by 
the Australian Geodetic Datum indicate that the accuracy of the 
geocentric orientation vector so determined from incomplete re- 
presentations of the Earth’s gravitational field, is dependent on the 
overall extent of the datum over which the comparisons are made. 

The lack of an astro-geodetic determination of the "geoid" due 
to an inadequate density of astro-geodetic stations over certain 
parts of the datum does not materially affect these determinations. 
The effect is only marginal if the station density of the area re- 
ferred to does not markedly differ from the average for the region. 
These tests indicate that sufficient astro-geodetic data exists for 
the connection of the North American Datum to a world geodetic 
system assembled from gravimetric considerations and referred 
to the geocenter with an accuracy of 0"2 or its linear equivalent 
in each component of the geocentric orientation vector. 


*On leave of absence from the University of New South Wales, Sydney, Australia. 
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THE GEOCENTRIC ORIENTATION VECTOR 
FROM LIMITED ASTRO-GEODETIC DATA 


1. INTRODUCTION 

The establishment of a world geodetic system based on gravimetric consid- 
erations has been mooted as a distinct possibility (Mather 1971b). The geodetic 
reference frame established in this manner has the considerable advantage of 
being directly related to the Earth's center of mass or geocenter , as a conse- 
quence of conditions implicit in the solution (ibid, p. 83). A number of technical 
problems which have to be coped with in the implementation of such a scheme are 
elaborated on by Mather (1972a). While the accuracy attainable from this type 
of determination is, in the short term, estimated as only being equivalent to that 
from the geometrical use of satellites using optical techniques, it is nevertheless 
of fundamental importance in the definition of the observation space as it provides 
a relationship between the geocenter and the observing stations at the surface of 
the Earth. The significance of this aspect is discussed in detail by Mather (1972b). 

The basic techniques can be described briefly as follows. The separation 
vector d is defined as the separation in Earth space between equivalent points P 
on the physically defined surface (e.g. , the Earth's surface; the geoid) and those 
on a system of reference used to relate the former, and defined with respect to 
the geocenter (e.g. , the telluroid ; the reference ellipsoid). Its component d 4 
along the axes Xj of a local Cartesian system of reference whose x 3 axis is coin- 
cident with the local normal to the reference surface, while the Xj and x 2 axes 
are oriented in the north and east directions respectively, defining the horizontal 
plane, are given by 


dj = R £; d 2 = R i?; d 3 = h s (1) 

where R is the mean radius of curvature of the equipotential surface of the ref- 
erence system passing through the point of reference. £ and?? are components 
of the deflections of the vertical in the meridian and prime vertical directions 
respectively, and h s is the evaluation above ellipsoid. 

h s is related to the height anomaly hj, the normal elevation h n , the ortho- 
metric elevation h 0 and the results of astro-geodetic levelling as shown in Fig- 
ure 1, by the following set of relations. 

P 

hs “ hd + ~ ^ 1 Ah 0 + Ng eo id = N + h 0 (2) 

geoid 
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Figure 1. Elevation Above Ellipsoid 


where N is the elevation of geoid above ellipsoid, the subscript geo id referring 
to the value at a point on the geoid from which the astro-geodetic levelling de- 
fining N a is commenced. The summation which appears after the second equality 
in equation 2, is taken over the route of the astro-geodetic levelling, h 0 being the 
increments in orthometric elevation. The results of astro-geodetic levelling, 
given as N a in equation 2, are defined by the relation 

P 

Na = - X) r d* (3) 

geoid 

where f is the component of the deflection of the vertical along the element of 
length di whose azimuth is a , it being assumed that the geoid has a linear gra- 
dient between the terminals of the element of length, f is defined by the relation 

J = £ cos a + tj sin a (4) 

Both h n and h 0 are linear interpretations of the difference in geopotential 
AW between the geoid and P which is an observed quantity. h n is defined on the 
reference system, being a displacement in space unoccupied by matter, in the 
gravitational field exterior to an ellipsoid of revolution whose bounding surface 
is an equipotential and with mass characteristics similar to that of the Earth 
(e.g. , Heiskanen & Moritz 1967, p. 171). h 0 on the other hand is a displace- 
ment along the local vertical in space occupied by the topography exterior to the 
geoid and within the Earth's surface. It follows that while Ah 0 is an observed 
quantity, h 0 is dependent on the model adopted for the Earth's topography ex- 
terior to the geoid. This has to be specified before h 0 and N have a unique defi- 
nition for locations in land areas. 
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Most of these distinctions are conveniently blurred in practice due to limita- 
tions on the accuracy of solutions as a consequence of the incomplete definition 
of the requisite data sets. Practical solutions are usually called representations 
of the "geoid" irrespective of whether they are merely the results of the astro- 
geodetic levelling operation which merely defines N a on its own, or if the deter- 
mination is gravimetric, obtained by the use of free air anomalies in Stokes' 
integral, given by 


N 


f 



m|a p I r rr 

- R — t~ + ^JJ f( ^ A s da 


(5) 


where Ag is the value of the free air anomaly representing the element of surface 
area da which is at an angular distance \ p from the point P at which Nf is deter- 
mined. f(\p) is Stokes' function whose formulation is well known (e.g. , ibid 1967, 
p. 94). W 0 is the potential of the geoid, U 0 that assigned to the surface of the 
reference ellipsoid and M j Ag [ the global mean of the free air anomaly. 

The free air geoid Nf so defined is a convenient quantity to work with as it 
embodies the total effect of observed gravity and further, Nf is the major con- 
tributor to both N and h d . While having significant magnitude which is approxi- 
mately 10% that of Nf, the correction terms could exhibit variations which can- 
not be ignored in certain instances. The expressions for these correction terms 
have been investigated and the required relations are 


N f = h d -N c - N-Nj 


( 6 ) 


N c is given by (Mather 1971b, p. 85) 


N P 


r2 rr 1 // dh \ N f 

= 27JJ7l( Rsin ' // d7 + (h P _h) )Tr- tan ^a 


da + o j fN c | 


if 


1 

2 



<f 


(7) 


r being the distance between the element da and the point of computation P. The 
terms to be integrated in equation 7 should be evaluated at da , h referring to the 
orthometric elevation, the subscript p indicating the value at P. tan |3, and tan 
/3 2 are the ground slopes in the north and east directions, while 


dh 

dr 


cos a tan j3j + sin a tan j3 2 


(8) 


3 



Also note that the use of the repeated greek subscript in equation 7 refers to sum- 
mation over two indices, while R and 7 in equations 5 and 7 refer to the mean 
radius and the mean value of gravity for the Earth model. For a detailed de- 
scription of Nj , see Mather 1971c, p. 84. The solutions described above are 
correct to the order of the flattening f and should suffice until such time as im- 
provements in the definition of the Earth's gravitational field warrant extensions 
of computation techniques. 

Fryer (1970, p. 155) has computed Nj using an approximate technique and 
his results indicate that its magnitude is generally small with maxima over the 
Himalayas and the Andes, but with rather slow attenuation. This calls for the 
exercise of caution when using the free air geoid alone for the representation of 
the gravimetric determination of the separation vector as has been done for an 
investigation covering the Australian Geodetic Datum (Mather 1971a). This is 
necessary as a significant distorting effect could be produced on the computed 
value of the geocentric orientation vector if any significant correlation were to 
exist between topographic variations and position. The region covered by the 
Indian Geodetic Datum would be an example of such an area. 

Such considerations would also apply, though to a lesser degre^, to the re- 
gion covered by the North American Datum (NAD) which has greater topographical 
variations than the Australian region. Consequently, it is possible to compare 
astro -geodetic solutions directly with the free air geoid over the latter datum 
without errors much in excess of ±1 m due to this cause and terrain variations 
produce insignificant effects on the computation of the geocentric orientation vec- 
tor (Mather, Barlow & Fryer 1971, p. 24). The NAD is therefore an excellent 
test region for the investigation of topographical effects in view of the quality of 
both the gravity and elevation data available in the area and the variety of solu- 
tions available for geocentric determinations from geometrical considerations, 
which would afford independent checks. 


2. THE MATCHING OF SOLUTIONS 

Solutions for the geocentric orientation vector can be visualized as the match- 
ing of imperfect determinations of the same surface using two different techni- 
ques. Such a description is somewhat simplistic in the light of the reasons given 
in the last section. The necessary equations are obtained by attributing the dif- 
ference in the values for the separation vector as obtained from gravimetry and 
astro-geodetic methods in turn, as being due to the latter being referred to an 
ellipsoid whose center is not at the geocenter, though its minor axis is parallel 
to that of the Earth's rotation for the epoch of observation. The solutions can 
be referred without difficulty to a reference figure of the same dimensions by 
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( 9 ) 


the standard relations (e.g. , Mather 1968, pp. 290 et seq) 


50 = df sin 20 + o {f 3 } ; 


5X = 0 (10) 

and 

5h s = da + a df sin 2 0 - v 50 tan 0 + o {f 3 a} (11) 

which link changes da in the equatorial radius a, df in the flattening f of the ref- 
erence ellipsoid with changes 50 in the latitude 0 , 5X in the longitude X , and 
6 h s in the ellipsidal elevation h s . v is the radius of curvature of the reference 
ellipsoid in the prime vertical section at the point considered. 

The geocentric orientation vector £, which is the space vector defining the 
displacement of the center of the regional ellipsoid with respect to the geocenter, 
as illustrated in Figure 2, can be represented by the components A£j in the 
Laplacian triad at the origin of the regional datum. The transformation between 
components dj a of the separation vector in the local Laplacian triad at any point 
of the regional datum, as determined from astro-geodetic considerations, can 
be related to those evaluated from gravimetry (dj g ) by the equations (Mather 1971a, 
p. 63 et seq) 


^ig ^ia + ft. ^ij < h o)j 


( 12 ) 


where the summation convention applies. Ajj are the elements of the array 


cos 0 O cos 0 sin 0 Q cos 0 

sin 0 sin AX 

+ sin 0 O sin 0 cos AX - cos 0 O sin 0 cos AX 


where 


A = 


- sin 0 Q sin AX cos AX 
sin 0 cos 0 


0 


sin 0 sin 0 


- cos 0 sin AX 


- sin 0 O cos 0 cos AX 


+ cos 0 O cos 0 cos AX 


(13) 


AX = X 0 -X, 
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Figure 2. The Geocentric Orientation Vector 


the subscript 0 referring to values at the origin of the regional datum at which 
the geocentric orientation parameters A£j are computed. The quantities hj are 
given by 

hj = - (p + h s ); h 2 = - (u + h s ); h 3 = 1, (14) 

where pis the radius of curvature of the reference ellipsoid in the meridian. 

Three types of solutions are possible for the quantities A£j. The first com- 
pares values of h s only which is equivalent to restricting the block of equations 
at 12 to those obtained when i = 3. The second is based on the comparison of 
values of £ and p which restricts the equations at 12 to those obtained when i = 1 
and i = 2. A third solution uses all three quantities £ , % j and h s . Expressions 
of gravimetric values of £ and r? which are directly comparable with astro- 
geodetic values are given in Mather 1971b, p. 87. 

Experiences based on tests carried out on the Australian Geodetic Datum in- 
dicate that solutions of type 1 and 3 give identical results. Solutions of type 2 
give slightly different results for A£j and A£ 2 (Mather 1971a, p. 72). These 
have been attributed to possible errors in the process used to deduce N a and 
hence h s from a density of astro-geodetic stations which is less than ideal. This 
conclusion is strengthened by the fact that the discrepancies between the values 
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of A£, and A£ 2 as determined by either one of types 1 and 3 or type 2, are 
smaller when the density of astro -geodetic stations is doubled. 

If the definition of the global gravity anomaly field is free from error, it 
follows that A£j can be obtained by a direct comparison of d g j and d a j at the or- 
igin alone. Such a technique cannot be resorted to in practice without signifi- 
cantly affecting the accuracy of the result, because the accuracy of the astro- 
geodetic determinations of 0 and X is restricted to ±0. 3 arcsec unless the direction 
of the vertical at the origin has been established by techniques more akin to those 
used at an observatory. A second difficulty which would be complicating is the 
effect any large scale prediction of gravity values has on the final results. 

Australia is a region where the gravity anomaly field was based on a station 
density which was adequate for geodetic calculations, affording representation on 
a 10 km grid with 80% coverage. Earlier tests indicated that serious errors are 
obtained when determining the geocentric orientation parameters A£, for an area 
this size if the values of A£ ; were obtained by a comparison of the components 
of the separation vector dj at locations which were restricted to a small part of 
the datum (Mather & Fryer 1970, p. 278). 

This is due to two causes. The first is the significant effect errors of pre- 
diction have on the final result. Such errors could be minimized in areas where 
no observed gravity is available over extensive regions by constraining predicted 
values to fit estimates based on those from satellite orbital analysis, as strengthened 
by combination with surface gravimetry (Mather 1971a, pp. 65 et seq). However 
this merely minimizes, rather than eliminates the error. 

The second is due to eri’ors in the representation of the largely unsurveyed 
ocean areas by low degree harmonics of the Earth's gravitational field obtained 
by the technique referred to in the pi'evious paragraph, as initially proposed by 
Kaula (1966). Errors in the higher harmonic coefficients which are determined 
with greater uncertainty, could seriously distort the geocentric orientation 
parameters A£j as determined from small areas, as illustrated in Figure 3. 

Tests on the orientation of the Australian Geodetic Datum (AGD) which covers 
an area equivalent, to that of the United States, indicated that the best results are 
obtained by comparing gravimetric and astro-geodetic values of the separation 
vector c[ at points which would provide the widest possible representation of the 
entire datum whose geocentric orientation vector is to be defined. The precision 
of such determinations is not dependent on the number of points at which com- 
parisons are made but rather on the extent of datum represented. Thus the use 
of 38 well chosen points for a determination of the AGD on this basis, gave es- 
sentially the same result as that obtained from comparisons at 693 points on an 
evenly spaced grid (Mather 1971a, p. 72). 
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Figure 3. Orientation Errors from Limited Fixes 



The value of the geocentric orientation vector obtained from such compari- 
sons, effectively allows for errors in the determination of individual values of d. 
It can be concluded that the errors in the final value of the geocentric orientation 
vector are a consequence of those harmonics in the global representation of the 
gravity field whose departures from linearity over the region covered by the 
datum are less than the errors inherent in the technique used for the 
determination. 


3. APPLICATIONS TO THE NORTH AMERICAN DATUM 

The North American Datum (NAD) is taken as the geodetic datum for the 
United States, Canada and Mexico (NAS Report 1971, p. 10). The surface area 
covered is approximately two and a half times that of the AGD. Determinations 
of the separation vector as defined in equation 1 are available at a widely dis- 
tributed network of points in the United States (Rice 1971) and southern Canada 
(Ney 1951; Corcoran 1967). While the density of the astro-geodetic stations over 
the rest of the datum is low, the U. S. Defense Mapping Agency's Topographic 
Command has produced a geoid map (i.e. , N a values) of North and Central 
America (Fischer et al 1967). Without detracting from the ability of the group 
working at this organization on the preparation of geoid maps from astro- 
geodetic data, which is second to none, it is always possible that the paucity of 
data can give either an over-smoothened result in these sparsely surveyed 
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regions or else a slight bias in the values of the geocentric orientation parameters 
A£j determined, in view of the dominating influence of d 3 on the solution, as men- 
tioned in Section 2. 

It is therefore preferable to provide the vital coverage for the large areas in 
northern Canada and Mexico by the deflections of the vertical alone. This is con- 
sidered to be necessary as the determination of the quantities A|j by restricting 
comparisons to the central region alone could bias the results for the reasons 
given in Section 2. 

The optimum technique for the determination of the geocentric orientation 
parameters A^ , and hence the geocentric orientation vector 0 is the following. 

Equation 13 is applied in two different ways depending on whether sufficient 
data were available for a reliable astro -geodetic determination of the ellip- 
soidal elevations h s (i.e., d 3 ). 

If it were, then three observation equations would be formed at each point, 
in terms of equation 13. Else, there would only be two, with the last row in 
the array on the right of equation 14 not being used. 

In this maimer, the final solution for the geocentric orientation parameters A^ 
will be representative of the comparisons over the entire datum. 

The use of this representation of the separation vector as determined from 
astro-geodesy was put to the test on the Australian Geodetic Datum where an 
astro -geodetic determination of h s is available for the entire region, with no sig- 
nificant (i.e. , < 01'2) bias between the deduced geo id solution and the observed 
deflections of the vertical (Mather et al 1971, p. 24). The results of a compari- 
son made over all the available astro-geodetic stations on the AGD is given in 
row 1 of Table 1. The Australian continental area extends in latitude between 
-10° N and -40° N, if Tasmania were excluded. 

It was decided to exclude the geoid solution from the determination of the 
geocentric orientation parameters A|j at latitudes above a certain minimum value 
0 m , and restrict the comparison of components of the separation vector at lati- 
tudes 0 (> 0 m ) to those of the deflections of the vertical only (i. e. , d 3 , d 2 ) . Com- 
parisons at points where 0 > 0 m would still be represented by the three observa- 
tion equations implied in equation 13. 

The adoption of the procedure outlined in the previous paragraph creates a 
situation akin to that existing at the present time for the NAD. The solutions 
obtained for the geocentric orientation parameters A£j as 0 m varied between 
-10° N and -34° N are set out in Table 1. 
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Table 1 


The Effect of Excluding the Geoid Solution 
from Determinations of the Geocentric Orientation Parameters 


0m 

No. of Astro -Geodetic 
Stations 

Geocentric Orientation Parameters 


0 > 0m** 

A £ ! (sec) 

A £ 2 (sec) 

A £ 3 (m) 

-10°N 

1084 

0 

-3.959 


8.30 

-13°N 

1071 

13 

-3.964 


8.31 

-16°N 

1022 

62 

-3.958 

-4.119 

8.29 

-19°N 

926 

158 

-3.961 

-4.113 

8.30 

-22°N 

750 

334 

-3.986 

-4.103 

8.38 

-25°N 

624 

460 

-3.988 

-4.112 

8.34 

-28°N 

444 

640 

-3.975 

-4.107 

8.24 

-31°N 

242 

842 

-4.012 

-4.088 

8.45 

-34° N 

85 

999 

-3.984 

-4.086 

7.92 


* Three observation equations per point 
**Two observation equations per point 


The resulting values of the orientation parameters show a stability which is 
certainly an order of magnitude better than the final accuracy attainable from the 
gravity data available at present. It should be pointed out that the density of 
astro -geodetic stations at latitudes greater than <fi m is no different to those at 
latitudes less than <p m as the astro-geodetic stations are evenly distributed over 
the Australian region. It may be reasoned that the results would therefore not 
be directly applicable to the North American situation. However, earlier tests 
on the AGD indicate that results of equivalent accuracy can be obtained from a 
few well distributed stations as is possible from a much higher density, though 
the errors in A£j are likely to be greater than those implied in the results from 
the Australian tests. This would still not bias the final values of the geocentric 
orientation parameters A£j in excess of the accuracy estimates given in the 
summary. 


4. CONCLUSION 

The relatively sparse distribution of astro -geodetic stations in the regions 
covered by the North American Datum, and lying outside the United States and 
southern Canada, need not be an inhibiting factor in defining the geocentric 
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orientation vector for the NAD with an accuracy equivalent to that attainable by 
the geometrical use of artificial Earth satellites based on optical techniques. A 
solution of adequate accuracy could be obtained if those regions not covered by 
networks of astro-geodetic levelling were represented by individual astro-geodetic 
stations at which the deflections of the vertical had been established, with the 
proviso that these stations afforded an evenly distributed representation of the 
regions. 

When considering the consequences of the relatively low density of astro- 
geodetic stations in these regions, on the values of the geocentric orientation 
parameters, it is instructive to recall that the Australian Geodetic Datum was 
established by setting the reference ellipsoid parallel to an estimate of the mean 
geoid slope across the region using only 150, but well spaced astro-geodetic sta- 
tions (Bomford 1967). The additional information provided by nearly eight times 
as many stations indicate that the consequent error between the desired condition 
and that achieved with this sparse but well distributed sample was only 0'.'2 in 
each of the orientation parameters at the origin (Mather et al 1971, p. 24). It 
is confidently estimated that the use of the technique proposed will introduce er- 
rors which are well within these figures. 
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